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Abstract 

The influence of cold nuclear matter on lepto-production of hadrons in semi- 
inclusive deep inelastic scattering is measured using the CLAS detector in 
Hall B at Jefferson Lab and a 5.014 GeV electron beam. We report the 
multiplicity ratios for targets of C, Fe, and Pb relative to deuterium as a 
function of the fractional virtual photon energy z transferred to the and 
the transverse momentum squared p\ of the . We find that the multiplicity 
ratios for are reduced in the nuclear medium at high z and low p^, with 
a trend for the transverse momentum to be broadened in the nucleus for 
large p\. 

Keywords: Hadronization, Hadron production, Deep inelastic scattering, 
Nuclei, Quarks 



Hadronization is the process through which partons, created in an ele- 
mentary reaction, turn into hadrons. An example for this process is lepto- 
production of hadrons in semi-inclusive deep-inelastic electron-nucleon scat- 
tering, where a parton is struck by a virtual photon; as the parton propagates 
away from the interaction point it picks up partner quarks (antiquarks) to 
form a hadron. The process of hadronization is changed by the nuclear 
medium. Nuclear modification of hadron production in deep-inelastic scat- 
tering (DIS) was first observed at SLAG [1] followed by EMC fl, E665 Q and 
more recently at HERMES 0-0] . The goal of these investigations is to study 
the mechanism of quark deconfinement followed by propagation and then re- 
confinement into a hadron, and in particular how the presence of cold nuclear 
matter affects the hadronization process. Although some rescattering of the 
quark/hadron in nuclear matter is expected, calculations 0, Isj for the HER- 
MES data suggest that the effects are more than just simple rescattering-the 
hadronization mechanism itself is changed due to the propagation through 
nuclear matter. In addition to understanding the fundamental features of 
hadronization, a detailed understanding of the nuclear modifications is im- 
portant for the interpretation of ultra-relativistic heavy ion collisions, as well 
as the study of neutrino oscillations using nuclear targets in the GeV energy 
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range. A recent review on the subject of hadronization and parton propaga- 
tion can be found in reference 

The observables reported in the present paper are the multiphcity ratios 
as a function of z = Eh/iy, the fraction of the virtual-photon energy car- 
ried by the hadron in the target rest frame, and its transverse momentum 
squared, p^. The transverse momentum is defined as the perpendicular com- 
ponent of the hadron's momentum measured with respect to the direction 
of momentum transfer in electron scattering. The multiplicity ratio R'X is 
defined as the ratio of the number of hadrons detected, normalized to the 
number of electrons measured in the DIS kinematics for a nuclear target with 
mass number A, divided by the same quantities for a deuterium target (D); 

where Nh is the number of hadrons and N^^^ is the number of scattered 
electrons detected in the specified kinematic bins. Here, is the negative 
of the four-momentum of the virtual photon squared and u is the energy of 
the virtual photon. For most of the explored kinematic range R'X shows a 
reduction from unity, although at low z and high it shows an enhancement. 
Hadron formation lengths, the characteristic distances over which hadrons 



form, can be extracted from [10 



In certain theoretical models, the transverse momentum is broadened for 
hadronization from nuclear targets as compared with the distribution mea- 
sured from deuterium, and the amount of broadening is sensitive to gluonic 
radiation by the quark before it evolves into a hadron jo], [ill, 12 1. By varying 



the nuclear radius, measurements of the broadening as a function of the kine- 
matic variables allow one to infer the length scale (known as the production 
length) over which the quark is deconfined. 

Since hadronization is a nonperturbative QCD process that as yet cannot 
be calculated from first principles, one must rely on models to interpret the 
observed nuclear dependence of hadronization. In particular, several hadron 
species are required to analyze the flavor dependence of the observables, 
to unravel the reaction mechanisms involved and to give insight into the 
fundamental principles governing hadron formation. Of special interest is 
the hadronization into strange mesons as compared with non-strange mesons. 
In the HERMES data the z dependence of R\ in the medium for 7r+, 
7[~ and 7r° is nearly identical for a given nucleus. However, R^ for is 
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systematically larger than for K~ . The difference for the charged kaons 
can be explained due to the larger nuclear absorption for K~ compared 
to . Also, note that pions have only non-strange valence quarks, 
has a non-strange valence quark while K~ has a non-strange antiquark and 
strange quark content. Thus, the production of a negative kaon requires 
further complicated string breaking mechanisms and fewer leading hadrons 
are expected to be K~ as compared to . 

While the HERMES measurements for the hadronization effects in 7r° pro- 
duction looked very similar to vr^, they did not report any hadronization 
results. Strangeness production can give a unique insight into the hadroniza- 
tion process, since the strange quark is likely to come directly from the 
"string-breaking" process Different hadron species, such as the K^^ are 
relevant to the study of the flavor dependence of hadronization and the mass 
dependence of quark energy loss. In this letter, we report for the first time 
on hadronization measured using the CLAS detector at Jefferson Lab. 
Our results show, like the HERMES results, that the hadronization mecha- 
nism is affected by the presence of the nuclear medium. Of course, the 
is a mixture of and , but just as hadronization dominates over 
K~ (presumably due to the non-strange valence quark in the K^) we expect 
that most of the signal is dominated by hadronization, with a smaller 
contribution from . 

The experiment was carried out using the CLAS detector fis'] at the 
Thomas Jefferson National Accelerator Facility. The data presented here 
are from the E02-104 [l^ experiment. An electron beam of 5.014 GeV was 
incident on targets of liquid deuterium, solid carbon, iron and lead. The 
targets were arranged in a mechanical structure that allowed one solid nuclear 
target to be in the beam at the same time as (and slightly downstream of) the 



liquid deuterium target pj|. The length of liquid target was about 2 cm and 
the separation distance between solid and liquid targets was approximately 4 
cm. This allowed for reduced systematic uncertainties in the ratios of hadrons 
measured by the CLAS detector. The scattered electron and the produced 
hadrons were detected in coincidence by the CLAS detector. Details of the 



experiment are described in Ref. [16 



The DIS kinematics covered by the present experiment are shown in 
Fig. [H where is plotted versus u and also versus Bjorken x = Q^/ (2Mi/), 
where M the nucleon mass. Constraints have been applied, limiting the 
kinematic variables to > 1.0 GeV^, > 4.0 GeV^ and y < 0.85, where 
is the squared invariant mass of the photon- nucleon system, y = u/E 



5 




1.5 2 2.5 3 3.5 4 4.5 5 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

V (GeV) X 

Figure 1: Distribution of versus v (left) and x (right), for DIS events from a 5.014 GeV 
electron beam on a liquid deuterium target, after the kinematic cuts described in the text. 

the energy fraction of the virtual photon in the target rest frame, and E is 
the beam energy. The constraint on y is apphed to reduce the magnitude 
of radiative corrections. The range in x depends on and extends from 
X = 0.1 at the lowest up to x = 0.55 at the highest . Though the sea 
quarks may have non-negligible contribution at the smallest x values probed 
in this experiment, for x greater than about 0.3 the virtual photon is likely 
to strike a valence quark. Note that the kinematics of the CLAS results are 
different from those of HERMES results, which had higher v and hence lower 
X, extending down to 0.023. 

The ifg 's were found from fits to the invariant mass distribution of tt+tt" 
pairs produced. After identifying the scattered electron, K^^ invariant mass 
spectra M(7r+, vr^) is formed from the combination of all charged pion pairs, 
detected in CLAS. For a given bin in and integrated over all electron 
kinematics and all other hadron variables, a clear peak in M(7r+, vr") is seen 
at the mass of the as shown in the top plots of Fig. |2]for both the liquid 
deuterium (left) and solid Pb (right) targets. Similarly, plots showing the 

peaks for a bin in after integrated over all electron kinematics and all 
other hadron variables are shown in the bottom plots of Fig. [2J Although 
substantial, the background from uncorrelated vr+vr^ pairs is smooth and 
is easily fit by a polynomial. The K^^ peak was fit using an unconstrained 
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Figure 2: Invariant mass of two oppositely charged pions, showing a peak at the mass, 
for the hquid deuterium target (left) and the solid Pb target (right) for selected kinematic 
bins: 0.3 < z < 0.4 (top) and 0.1 < < 0.2 GeV^ (bottom). The x-bin size is 2 MeV. 
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Gaussian. The statistical uncertainties are determined directly from the fit 
parameters. 

Results for R'\ as a function of z and of p"^ for hadronization are shown 
in Figs. |3] and HJ respectively, for the carbon, iron and lead targets. The 
kinematic dependence on a given variable is extracted while integrating over 
all other variables within the acceptance of the experiment. Results for vs. 
z for the charged kaons from the HERMES experiment are also shown in left 
panel of Fig. [31 For comparison, we have plotted only the results for the xenon 
target from the HERMES experiment; results for several nuclei and different 
hadrons are available in Ref. For the p'^ dependence studies, the data 
presented are limited to 0.3 < z < 0.8 to reduce the possible contamination 
from target fragmentation and from exclusive channels. The ratios of 
have been normalized by the number of DIS electrons for each target and 
corrected for geometrical detector acceptance for each target as calculated 
by Monte Carlo simulations using a PYTHIA event generator ly, |r7| . The 
acceptance correction was applied on a bin-by-bin basis. No corrections have 
been applied to the data for QED radiative effects, which were too small to 
measure with the present statistics. However, our studies showed that these 
corrections are less than 5% for the multiplicity ratios and this is included 
in the systematic uncertainty. 

Several sources contribute to the overall systematic uncertainty in the 
multiplicity ratios. Details of the systematic uncertainty studies are de- 



scribed in Ref. [18j. The systematic uncertainty in the acceptance correction 
was found to be less than 5% for the z dependence (except for the lowest z 
bin where this uncertainty is estimated to be 16%) and less than 8% for the 
dependence. The systematic uncertainty associated with particle identi- 
fication for the multiplicity ratios was found to be less than 3%. To account 
for the possible systematic uncertainty in the fitting procedures we have as- 
signed a 4% uncertainty in R\ vs. z results and a 5% uncertainty in R\ vs. 
p'^. The total uncertainty is obtained by adding systematic and statistical 
uncertainties in quadrature, and are shown by the longer error bars in Figs. [3] 
andm 

Fig. [3] shows that, for z > 0.3, fewer are seen from nuclear targets 
(normalized by the number of DIS events) than from deuterium while at z < 
0.3 the ratio shows a trend to be bigger than unity. At lower z, effects such 
as nuclear rescattering, including (vr, K) reactions, and target fragmentation 
are likely to contribute to the number of K'^ measured. The same effect at 
low z for hadronization has already been noted by HERMES fl. The 
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Figure 3: Left panel shows the CLAS data from the EG2 experiment for the hadronic 
multiplicity ratio for versus z along with HERMES results [G] for charged kaons for 
the Xe nuclei. Right panel shows the CLAS results along with the results derived from 
calculations of Ref. [3] for Pb (solid line), Fe (dashed line) and C (dotted line). Data 
points are shown with different symbols as given in the legend. The inner error bars 
represent the statistical uncertainty, while the outer ones show the total uncertainty. For 
clarity, the HERMES points are shifted by +0.02 units in z and the C and Pb points are 
shifted by -0.01 and +0.01 units in z, respectively. 



trend of the CLAS hadronization data is similar to that of the HERMES 
results, where more attenuation is seen as z increases, with larger attenuation 
factors for the heaviest target. The present uncertainties are large due to the 
limited statistics, however the data agree reasonably well with predictions [3] 
as shown in right panel of Fig. [31 

The theoretical calculations shown in Fig. [3] and presented in Ref. j3] 
are carried out in the framework of a probabilistic coupled-channel transport 
model based on the Boltzmann-Uehling-Uhlenbeck equation, which allows 
for detailed treatment of the final-state interactions beyond simple absorp- 
tion mechanisms. It starts with the Lund string fragmentation picture as 



embodied by PYTHIA [17| and extracts two points for each event: (i) the 



location in space where the struck quark breaks the string, and (ii) the loca- 
tion where the quark joins with an antiquark to form a meson. This is done 
event-by-event, in the context of a Monte Carlo calculation. Medium effects 
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are included by assuming a phenomenological transition from a pre-hadronic 
scattering cross section to a full cross section constrained by hadronic beam 
data. A function linear in the formation and production times is found to 
best fit the HERMES and EMC data (see Equation 3 and Fig. 1 of Ref. 0). 
Note that the calculations presented in [3] are for K'^ and For a given 
nucleus, we show the average of the above calculations to compare with the 
K^^ data. 

In some theoretical models joj , the decrease in R\ with z can be explained 
by the struck parton losing energy before picking up a partner antiquark, on 
the way to forming a hadron. One mechanism of energy loss is gluon radiation 



in the nuclear medium, although other mechanisms such as absorption [19 
are also possible. This energy loss reduces the hadron energy and hence 
reduces z for hadronization from a nucleus. Because the fragmentation func- 
tions have a steep dependence on z, the ratio R'X generates a steady decrease 
with z. Other models assume a decreased formation length for the struck 
quark to pick up a partner antiquark, followed by rescattering and absorption 
of the resulting hadron. 

The results for at CLAS kinematics show for z > 0.6 with the Pb 
target to be below the HERMES results for charged kaons from xenon. One 
must be careful when comparing these two data sets because of the different 
kinematic ranges; however, the CLAS kinematics provide an additional point 
of constraint for theoretical models of hadronization. Note that neutral kaons 
have a smaller reaction cross section than pions, and hence rescattering effects 
are minimized for hadronization. In particular, the data should be less 
affected by hadronic final-state interactions than pions, and more sensitive 
to QCD effects (such as gluon radiation) as the quark propagates through 
cold nuclear matter. 

Fig. m shows the multiplicity ratios as a function of the transverse mo- 
mentum squared. At high p^, there are very few K'^ events, leading to 
large statistical uncertainties. However, at small p"^ the statistics are rea- 
sonable and there is a clear target-dependence to the nuclear attenuation for 
p'^ < 0.1 GeV^. The multiplicity ratios then increase as p'^ increases and 
exceed unity at around p^=0.45 GeV^. This crossover is at lower than 
for the HERMES results. Furthermore, the CLAS data have a steeper slope 
with than those measured by HERMES. Both CLAS and HERMES data 
show that R\ levels off for lighter targets at low yet continues to decrease 
for heavier targets at low p"^. 

The increase in R'X with p'^ was seen previously by EMC [2I] and HER- 
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Figure 4: CLAS data for the hadronic multiplicity ratio for as a function of for 
0.3 < z < 0.8. Error bars are as in Fig. [3] For clarity, the C and Pb points are shifted by 
-0.01 and +0.01 units in p^, respectively. 



MES [sl, |6[. This is known as the Cronin effect 20], in which rescattering of 
either the parton or the hadron, as well as hadronic final state interactions, 
pushes events up to larger p'^. These effects occur preferentially in nuclear 
targets. The HERMES data for pions [gI, Esf shows that the Cronin effect 
diminishes for high- 2; events, which is consistent with predictions from parton 
rescattering [s!]. The question whether this also happens for hadroniza- 
tion must await better statistics, which will be available after the CLAS 
detector is upgraded in a few years for 12 GeV beam energies. 

In summary, we have extracted for the first time the multiplicity ratios 
for using semi- inclusive deep inelastic scattering with detection of the 
Kg using the CLAS detector with a 5.014 GeV electron beam incident on 
both deuterium and nuclear targets. Many systematic effects cancel in the 
ratio R\, in part because both deuterium and nuclear targets were exposed 
to beam simultaneously and were separated by only a few centimeters at the 
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center of CLAS. The z and dependencies of B!\ follow the trends seen in 
previous data by HERMES Collaboration at higher v and Q^, suggesting that 
the general mechanism of hadron formation are similar at CLAS kinematics. 

In order to fully understand the fundamental principles of the hadron for- 
mation as well as the flavor dependence of the observables, a wide spectrum 
of both baryons and mesons are required. Some of these channels, including 
more precise pion data, are currently under investigation by the CLAS col- 
laboration. Ultimately, the interpretation of these data will require further 
theoretical development. 
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